Spin susceptibility and magnetic short— range order in the Hubbard model 
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The uniform static spin susceptibility in the paraphase of the one-band Hubbard model is 
calculated within a theory of magnetic short-range order (SRO) which extends the four-field slave- 
boson functional-integral approach by the transformation to an effective Ising model and the self- 
consistent incorporation of SRO at the saddle point. This theory describes a transition from the 
| paraphase without SRO for hole dopings 8 > 8 C2 to a paraphase with antiferromagnetic SRO for 

8 C1 < 5 < 8c 2 . In this region the susceptibility consists of interrelated 'itinerant' and 'local' parts 
and increases upon doping. The zero-temperature susceptibility exhibits a cusp at 8 C2 and reduces 
to the usual slave-boson result for larger dopings. Using the realistic value of the on-site Coulomb 
repulsion U — 8t for La,2-s Sr{Cu04, the peak position (<5 C2 = 0.26) as well as the doping dependence 
reasonably agree with low-temperature susceptibility experiments showing a maximum at a hole 
doping of about 25%. 
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Among the most striking features of high-T c superconductors in the normal state, the unconventional magnetic 
properties have attracted increasing attention As revealed by neutron scattering M and nuclear magnetic reso- 
nance H experiments, in the metallic state there exist pronounced antiferromagnetic (AFM) spin correlations which 
are ascribed to strong Coulomb correlations within the Cu02 planes. Knight-shift B and bulk measurements [(i]-|6) of 
the spin susceptibility x(T,S) in La2- < 5Sr,5Cu04 show a maximum in the doping dependence as well as, for moderate 
hole doping (<5 < 0.21), in the temperature dependence, where the temperature of the maximum decreases with in- 
creasing doping. Such a behaviour, also observed in YBa2Cu306+ y (y < 0.92) f7)||, may be qualitatively understood 
I i as an effect of AFM short-range order (SRO) which decreases with increasing doping and temperature. 

Up to now there have been only a few attempts, based on one-band [PHf^] and three-band |u| correlation 
models, to describe the unusual doping and temperature dependence of the normal-state susceptibility. In the t-t'-J 
model, a maximum in x was obtained for the Pauli susceptibility of a strongly renormalized quasiparticle band jj| 
or for the RPA slave-boson susceptibility jf(J showing a cusp in the temperature dependence at the transition to 
the singlet RVB state. In the one-band Hubbard model, a maximum in the doping dependence of \ was found by a 
semi-phenomenological weak-coupling approach JTlJ or by the composite operator method jjjfl . The role played by 
SRO in explaining the normal-state susceptibility was investigated on the basis of the three-band Hubbard model 113] 
by means of a slave-boson CPA theory which, however, is self-consistent only at the single-site level and does not 
hold at very low temperatures. To improve the treatment of SRO in the paraphase being valid also at T — 0, in a 
previous communication j^J, hereafter referred to as I, we have presented the main features of a theory of magnetic 
SRO in the one-band Hubbard model based on the scalar four-field slave-boson (SB) approach |T5| ]. In I we have 
focused on the stability of magnetic long-range order (LRO) versus SRO, where magnetic LRO phases are found to 
make way to a paraphase with SRO in a wide doping region. 

In this Brief Report we extend our theory by the inclusion of an external magnetic field h and by the calculation 
of the uniform static spin susceptibility x i n the paraphase, where special care is taken to the influence of SRO. 

Following the lines indicated in I, the action of the SB functional integral for the partition function of the 2D 
Hubbard model is expressed in terms of the SB fields mi, £i, rij, i^, di and d* Jf6|. To treat the fluctuations of the 
local magnetizations m t and the internal magnetic fields £j we write m i = m,-^, ^ = £ i s i (s i — ±) and make the 
ansatz — > b Si for the magnetic amplitudes b € {fh, £} and the charge degrees of freedom b 6 {n, v, d = d*}. We 
transform the free-energy functional W to an effective Ising model in the nearest-neighbour pair ((ij)) approximation 
and obtain 
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The single-site and two-site fluctuation contributions <P aa = ( ! > ia {a. i )\ a =a and ^aa'a — ^(ij)cr( a i, a j) \ <\== > respec- 
tively, are given by 

<*V = i y dw/(w - M)Im In [1 - G&^fa)] , (5) 

# fe> = i y dtjf{w- n)1m In 1 - G {t:j)a T ja (a J )G {:jl)a T l<J (a i ) . (6) 

In (^), G°j a (uj) is the uniform paramagnetic (PM) Green propagator, and the scattering matrix Tj CT = 
V ia (1 — G° ia V ia ) ~ 1 is expressed in terms of the local perturbation 

V ia (au «) = 7Z~ ^2" { IK,) 2 - (4) 2 ] [w - v° + <x(f + ft)] + (^) 2 - ^° - a(a^ aj - ?)] } , (7) 
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and the superscript "o" refers to PM saddle-point values. By the functional (Q) we determine the saddle point for 
all Bose fields b a = {fh a ,£ a , n a , v a , d a } in the external field h, where, in the spirit of I, the SRO is self-consistently 
incorporated within the Bethe cluster approximation (taking into account only the nearest-neighbour SRO). As found 
in I, in the h = limit (fh a = m) the self-consistent calculation of the effective Ising-exchange integral J as function 
of the interaction strength U and the hole doping S — 1 — n yields two possible paraphases ((s,) = 0): (i) the paraphase 
without SRO (PM; J = 0, m = 0) and (ii) the paraphase with antiferromagnetic SRO (SRO-PM; J < 0, fa > 0). 
The uniform static spin susceptibility x(T, 5) has to be calculated according to 
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where m a — fh a a, W a = W a (h, h* , J) is the probability for the Ising spin a at the central site of the Bethe cluster, 
and h* is the effective Bethe field. The first term in (Jo[) describes the change of the magnetization-amplitude with the 
applied magnetic field and gives mainly the 'itinerant' contribution to \- The second term describes directional fluctu- 
ations of the local magnetizations and is called the 'local' contribution being finite only in the SRO-PM phase. Note 
that the 'itinerant' and 'local' properties are interrelated and determine both contributions to the spin susceptibility. 
In the PM and SRO-PM phases we have calculated the doping dependence of the zero-temperature susceptibility in 
the 2D Hubbard model (being finite in contrast to the theory of Ref. [jl3)) in a completely self-consistent way, where 
in the tedious numerical evaluation of the integrals (S) and (g) and of their derivatives particular attention has to be 
paid to the analytical behaviour of the complex logarithm. 
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FIG. 1. Uniform static spin susceptibility as a function of doping at T = 0. The theoretical result obtained for the 2D 
Hubbard model at U/t = 8 and t = 0.3 eV (solid) is compared with the spin contribution (x) to the (corrected) experimental 
susceptibility on Laa-«Sr«Cu04 at T = 50 K 

Figure 1 shows our result without any fit procedure using the commonly accepted value U/t — 8 for the Hubbard 
model applied to high-T c cuprates |17| . As stated in I, in the region 6 < U/t < 12, there occurs a first-order 
(l,l)-spiral^SRO-PM transition at 5 C1 and a SRO-PM^PM transition of second order at S C2 . In the PM phase 
(6 > 6 C2 ) the SB band-renormalized Pauli susceptibility has a pronounced doping dependence in two dimensions 
and agrees with the static and uniform limit of the dynamic spin susceptibility derived, within the spin-rotation- 
invariant SB scheme |D| , from the Gaussian fluctuation matrix at the PM saddle point |n| . In the SRO-PM phase 
(<5 Cl < 5 < 5 C2 ), the Pauli susceptibility is suppressed due to the SRO-induced spin stiffness against the orientation of 
the local magnetizations along the homogeneous external field. Accordingly, at S C2 a cusp in x(0, S) appears. Since, 
for 5 Cl < 6 < 5 C2 , \J\ decreases with increasing S fl4|| , the susceptibility increases upon doping. 

The peak in x(0,<5) only appears at sufficiently high ratios U/t > 6, for which a SRO-PM?^PM transition may 
occur. According to the phase diagram, given in Fig. 2 of I, in the region 6 < U/t < 12 the SRO-PM^PM transition 
shifts to higher doping values with increasing U /t. Correspondingly, the peak position in x(0, S) reveals the same U/t 
dependence. 

In Fig. 1 we have also depicted the spin contribution to the magnetic susceptibility of La2-^Sr5Cu04 at 50 K 
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obtained from the experimental data on the total susceptibility [|j by subtracting the diamagnetic core (—9.9 x 
10~ 5 emu/mol) and Van Vleck (2.4 x 1CP 5 emu/mol) contributions which, according to Ref. ||, can be taken to 
be independent of doping and temperature over the limited parameter region studied here. As Fig. 1 shows, the 
experimentally observed pronounced maximum at a hole doping of about 25% is reproduced very well by our theory 
yielding the peak position at d C2 — 0.26 (U/t = 8). Moreover, the qualitative doping dependence of x reasonably 
agrees with experiments. Of course, it could not be expected that our approach based on the simple (single-band) 
Hubbard model yields the correct magnitude of x f° r I^-iSraCuO^ Especially, concerning the low-doping limit 
S — > 5 = 0.04, the theoretical susceptibility is much too low as compared with experiments. This deficiency 
may be explained as follows. For 6 — and large U/t values, the Hubbard model is equivalent to the Heisenberg 
antiferromagnet with the exchange interaction J = 4t 2 /U. In this model, the spin susceptibility at T = has a 
finite value proportional to J -1 [ p0[ which is due to the existence of transverse spin fluctuations. However, our scalar 
four-field SB approach to the spin susceptibility in the presence of SRO implies the transformation of the free-energy 
functional to an effective Ising model describing longitudinal fluctuations only. Since the 'local' contribution to x is of 
Ising-type, we get a too small susceptibility in the low-doping limit which, however, is finite due to the interrelation 
to the 'itinerant' contribution to x- Therefore, we suggest that a theory of SRO based on the spin-rotation-invariant 
SB scheme fTsj ] and resulting in an effective Heisenberg-model functional may improve the results in the magnitude 
of Xi m particular at low doping levels. 

Finally, we notice that the increase of the susceptibility upon doping obtained within our theory for moderate 
Coulomb repulsions (U/t > 6) is in qualitative accord with recent QMC data and with the approaches of Refs. 

and p^ |. However, in those works a maximum in the spin susceptibility was found even at a smaller coupling 
W/t = 4). 

From our results we conclude that the concept of magnetic SRO in strong-correlation models may play the key 
role in the explanation of many unconventional properties of high-T c compounds. The theory may be extended in 
several directions. As discussed above a spin-rotation-invariant theory of SRO may improve the agreement of the spin 
susceptibility with experiments. Furthermore, as motivated by neutron scattering experiments M probing the AFM 
correlation length over several lattice spacings, the effects of a longer than nearest-neighbour ranged SRO (which may 
be described beyond the nearest-neighbour pair approximation) should be investigated. 

This work was performed under the auspices of Deutsche Forschungsgemeinschaft under project SF-HTSL-SRO. 
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